Search Results

Detailed chemistry represents a fundamental pre-requisite for a realistic simulation of combustion process in Diesel engines to properly reproduce ignition delay and flame structure (lift-off and soot precursors) in a wide range of operating conditions. In this work, the authors developed reduced mechanisms for n-dodecane starting from the comprehensive kinetic mechanism developed at Politecnico di Milano, well validated and tested in a wide range of operating conditions [1]. An algorithm combining Sensitivity and Flux Analysis was employed for the present skeletal reduction. The size of the mechanisms can be limited to less than 100 species and incorporates the most important details of low-temperature kinetics for a proper prediction of the ignition delay. Furthermore, the high-temperature chemistry is also properly described both in terms of reactivity and species formation, including unsaturated compounds such as acetylene, whose concentration controls soot formation.

The automotive industry is striving to adopt model-based engine design and optimization procedures to reduce development time and costs. In this scenario, first-principles gas dynamic models predicting the mass, energy and momentum transport in the engine air path system with high accuracy and low computation effort are extremely important today for performance prediction, optimization and cylinder charge estimation and control. This paper presents a comparative study of two different modeling approaches to predict the one-dimensional unsteady compressible flow in the engine air path system. The first approach is based on a quasi-3D finite volume method, which relies on a geometrical reconstruction of the calculation domain using networks of zero-dimensional elements. The second approach is based on a model-order reduction procedure that projects the nonlinear hyperbolic partial differential equations describing the 1D unsteady flow in engine manifolds onto a predefined basis.

A correct prediction of the initial stages of the combustion process in SI engines is of great importance to understand how local flow conditions, fuel properties, mixture stratification and ignition affect the in-cylinder pressure development and pollutant formation. However, flame kernel growth is governed by many interacting processes including energy transfer from the electrical circuit to the gas phase, interaction between the plasma channel and the flow field, transition between different combustion regimes and gas expansion at very high temperatures. In this work, the authors intend to present a comprehensive, multi-dimensional model that can be used to predict the initial combustion stages in SI engines. In particular, the spark channel is represented by a set of Lagrangian particles where each one of them acts as a single flame kernel.

In the most elementary treatment of plane-wave reflection at the open end of a duct system, it is often assumed that the ends are pressure nodes. This implies that pressure is assumed as a constant at the open end termination and that steady flow boundary condition is supposed as instantaneously established. While this simplifying assumption seems reasonable, it does not consider any radiation of acoustic energy from the duct into the surrounding free space; hence, an error in the estimation of the effects of the flow on the acoustical response of an open-end duct occurs. If radiation is accounted, a complicated three-dimensional wave pattern near the duct end is established, which tends to readjust the exit pressure to its steady-flow level. This adjustment process is continually modified by further incident waves, so that the effective instantaneous boundary conditions which determine the reflected waves depend on the flow history.

In the paper several application techniques of MonteCarlo (MC) method applied to thermal analysis of space vehicles are presented. Although these methods are widely used in other engineering domains, their introduction to the thermal one is quite recent and not fully developed in the industrial practice. This paper aims at showing that, even without demanding computation resources (all what presented has been obtained with a single processor PC) MonteCarlo analysis techniques, in a preliminary design phase, can support and integrate engineering judgment of the thermal designer. In particular, it is exploited the applicability of the method to reduced thermal models, with a clear advantage in terms of computation time. An original approach is proposed, and results are shown. The papers shows the applicability of the MC method to the case when experimental data of the uncertain parameters are available, using the bootstrap re-sampling techniques.

The current development to set up an automatic procedure for automatic mesh generation and automatic mesh motion for internal combustion engine simulation in OpenFOAM®-2.2.x is here described. In order to automatically generate high-quality meshes of cylinder geometries, some technical issues need to be addressed: 1) automatic mesh generation should be able to control anisotropy and directionality of the grid; 2) during piston and valve motion, cells and faces must be introduced and removed without varying the overall area and volume of the cells, to avoid conservation errors. In particular, interpolation between discrete fields is frequent in computational physics: the use of adaptive and non-conformal meshes necessitates the interpolation of fields between different mesh regions. Interpolation problems also arise in areas such as model coupling, model initialization and visualisation.

The s.i. combustion process and its corresponding pollutant formation are investigated by means of a quasiD approach and a CFD model. This work has been motivated by the need to better understand the reliability of such models and to assess their accuracies with respect to the prediction of engine performances and emissions. An extended dissertation about the fundamental mechanisms governing the pollutant formation in the turbulent premixed combustion which characterizes the s.i. engines is given. The conclusion of such analysis is the definition of a new reduced chemical scheme, based on the application of partial-equilibrium and steady-state assumptions for the radicals and the solution of a transport equation for each specie which is kinetically controlled. For this purpose the CFD code OpenFOAM [1, 2, 3] and the thermo-fluid dynamic code GASDYN [4, 5] have been applied and enhanced.

In the last years the number of electronic controllers of vehicle dynamics applied to chassis components has increased dramatically. They use lookup table of the primary order vehicle global parameters as yaw rate, lateral acceleration, steering angle, car velocity, that define the ideal behavior of the vehicle. They are usually based on PID controllers which compare the actual behavior of every measured real vehicle data to the desired behavior, from look up table. The controller attempts to keep the measured quantities the same as the tabled quantities by using ESP, TC (brakes and throttle), CDC (control shocks absorbers), EDIFF(active differential) and 4WS (rear wheels active toe). The performances of these controls are good but not perfect. The improvement can be achieved by replacement of the lookup tables with a fast vehicle model running in parallel to the real vehicle.

The air fuel mixing process of a small direct injection (d.i.) diesel engine, equipped with two different re-entrant combustion chambers and two nozzles having unlike spray angles, has been studied by integrated use of in-cylinder laser Doppler velocimetry (LDV) measurements, engine tests, and KIVA simulations. The LDV measurements have been carried out in an engine with optical access motored at 2200 rpm. The engine tests have been performed on a similar engine at the same speed, at fixed start of combustion, and different air-fuel ratio. The KIVA-II simulations have been made using as initial conditions the parameters determined by LDV and engine tests. The re-entrant bowl with higher levels of air velocity and turbulent kinetic energy at the time of injection gives the best performance. The nozzle having a spray angle of 150° which injects the fuel into the regions at higher turbulent kinetic energy lowers the smoke emission levels.

The influence of spray-wall interaction on air entrainment in an unsteady non-evaporating diesel spray was studied using laser Doppler anemometry. The spray was injected into confined quiescent air at ambient pressure and temperature and made to impact on a flat wall. The air velocity component normal to a cylindrical surface surrounding the spray was measured during the entire injection period, allowing to evaluate the time history of the entrained air mass flow rate. The influence of wall distance and spray impingement angle on air entrainment characteristics has been investigated and the results indicate that the presence of a wall increases the entrained mass flow rate in the region close to the surface, during the main injection period. Normal impingement appears to produce stronger effects than oblique incidence at 30 and 45 deg. A qualitative explanation of the results is also proposed, based on the drop-gas momentum exchange mechanism.

The paper illustrates an experimental and numerical investigation of the flow generated by an intake port model for a heavy duty direct injection (HDDI) Diesel engine. Tests were carried out on a steady state air flow test rig to evaluate the global fluid-dynamic efficiency of the intake system, made by a swirled and a directed port, in terms of mass flow rate, flow coefficients and swirl number. In addition, because the global coefficients are not able to give flow details, the Laser Doppler Anemometry (LDA) technique was applied to obtain the local distribution of the air velocity within a test cylinder. The steady state air flow rig, made by a blower and the intake port model mounted on a plexiglas cylinder with optical accesses, was assembled to supply the actual intake flow rate of the engine, setting the pressure drop across the intake ports atûP=300 and 500 mm of H2O.

In-cylinder cycle-resolved LDV measurements have been made in a diesel engine having a high-squish re-entrant combustion chamber with compression ratio of 21:1. The engine has been motored in the range of 1000 to 3000 rpm thanks to the use of self-lubricating seeding particles. Conventional ensemble-averaging and filtering techniques have been used for analyzing instantaneous velocity data obtained at two points along a diameter located in a horizontal plane at 5 mm below the engine head. The dependence of the mean motion and turbulence on engine speed has been evaluated. The effect of cut-off frequency selection on turbulence values has been also analyzed. Moreover, the Kolmogorov's -5/3 power domain has been investigated in detail by spectral analysis on the instantaneous velocity data.

In-cylinder measurements of turbulent integral length scales, carried out during the last 60 degrees of the compression stroke at 600 and 1,000 rpm by a two-probe volume LDV system, were used to assess the capability of the k-ε model used in KIVA-II code. The objective of the paper is to address the following question: what is the most reasonable definition of turbulent length scale in the k-ε model for engine applications? The answer derived from the comparison between KIVA predictions and experiments that showed a fair agreement between the computed turbulent length scale and the measured lateral integral length scale. The agreement is a result of proper choice of the initial swirl ratio and turbulent kinetic energy at inlet valve closure (IVC) by taking into account the LDV measurements and the value of the constant Cμε in the k-ε model equations that relates the turbulent length scale to k and ε.

Tangential component of velocity and turbulence were measured in three locations in the re-entrant combustion chamber of a motored single-cylinder d.i. Diesel engine (0.435 liter, 21:1 compression ratio) using a Laser Doppler Velocimetry system. Moreover, a modified LDV system with two-probe volume was used to measure directly lateral integral length scales of the velocity tangential component at two engine speeds. The measurements were made on a horizontal plane at 5 mm below the engine head from 100 degrees before TDC to 60 degrees after TDC of both the compression and expansion strokes. The engine was motored at 1,000 and 1,500 rpm respectively. An ensemble-averaging technique was performed to analyze the instantaneous velocity information supplied by two Burst Spectrum Analyzers. The lateral integral length scale was obtained from the integral of the spatial correlation coefficient of the velocity fluctuation for different separation.

Computations were carried out to simulate in-cylinder flow field and mixture preparation of a small port scavenged direct-injection two-stroke spark-ignition engine using a modified version of KIVA-3 code. Simulations of the interaction between air flow and fuel were performed on a commercial Piaggio (125 cc) motorcycle engine modified to operate with a hollow-cone injector located in different positions of the dome-shaped combustion chamber. The engine has a large exhaust port and five smaller transfer ports connecting the cylinder to the crankcase. The numerical grid of this complex geometry was obtained using an IBM grid generator based on the output of engine design by CATIA solution. To take into account the rapid distortion of flow, the standard k-ε turbulence model in KIVA-3 was replaced by the RNG k-ε model.

The present paper aims at developing a general method to estimate integral and microtime scales of turbulent in-cylinder flow field in reciprocating engines. The ensemble average technique was used to compute the integral time scale from the single point time autocorrelation function, whereas the microtime scale, representative of the most rapid changes that occur in the fluctuation, was computed as the intercept of the parabola that matches the autocorrelation function at the origin. Further, the microtime scale was also estimated by spectral analysis through the energy spectral density function of the ensemble turbulent fluctuation and the results obtained by the two methods were compared. The procedures were applied to the tangential component of the instantaneous velocity data collected, at different engine speeds (1,000, 1,500, 2,000 rpm), within a motored d.i. diesel engine equipped with a re-entrant combustion chamber, using the Laser Doppler Anemometry (LDA) technique.

A comprehensive model for sprays emerging from high pressure swirl injectors for GDI engine application has been developed. The primary and secondary atomization mechanism as well as the evaporation process both in standard and superheated conditions are taken into account. The spray modelling after the injection is based on the Liquid Instability Sheet Atomization (LISA) approach, modified to correctly predict the liquid sheet thickness at the breakup length. The effect of different values of the superheat degree on evaporation and impact on the spray distribution and fuel-air mixing is analyzed. Comparisons with experimental data show good agreements under atmospheric conditions and with different superheated degrees, while some discrepancies occur under higher ambient pressures.

Aim of this work is to present and discuss the possibility and the limits of two zone models for spark-ignition engines using a detailed kinetic scheme for the characterization of the evolution of the air-fuel mixture, while an equilibrium approach is used for the burnt zone. Simple experimental measurements of knocking tendency of different fuels in ideal reactors, such as rapid compression machines and shock tube reactors, cannot be directly used for the analysis of octane numbers and sensitivity of hydrocarbon mixtures. Thus a careful investigation is very useful, not only of the combustion chamber behavior, including the modelling of the turbulent flame front propagation, but also of the fluid dynamic behavior of the intake and exhaust system, accounting for the volumetric efficiency of the engine.

One of the most important phases in the development of direct-injected diesel engines is the optimization of the fuel spray evolution within the combustion chamber, since it strongly influences both the engine performance and the pollutant emissions. Aim of the present paper is to provide information about mixture formation within the combustion chamber of a heavy-duty direct injection (HDDI) diesel engine for marine applications. Spray evolution, in terms of tip penetration, is at first investigated under quiescent conditions, both experimentally and numerically, injecting the fuel in a vessel under ambient temperature and controlled gas back-pressure. Results of penetration and images of the spray from the optically accessible high-pressure vessel are used to investigate the capabilities of some state-of-the-art spray models within the STAR-CD software in correctly capturing spray shape and propagation.

The continuous pursuit of higher combustion efficiencies, as well as the possible usage of synthetic fuels with different properties than fossil-ones, require reliable and low-cost numerical approaches to support and speed-up engines industrial design. In this context, SI engines operated with homogeneous ultra-lean mixtures both characterized by a classical ignition configuration or equipped with an active prechamber represent the most promising solutions. In this work, for the classical ignition arrangement, a 3DCFD strategy to model the impact of the ignition system type on the CCV is developed using the RANS approach for turbulence modelling. The spark-discharge is modelled through a set of Lagrangian particles, whose velocity is modified with a zero-divergence perturbation at each discharge event, then evolved according to the Simplified Langevin Model (SLM) to simulate stochastic interactions with the surrounding gas flow.